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A recent approach to arrival times used the fluorescence of an atom entering a laser illuminated 
region and the resulting arrival-time distribution was close to the axiomatic distribution of Kijowski, 
but not exactly equal, neither in limiting cases nor after compensation of reflection losses by normal- 
ization on the level of expectation values. In this paper we employ a normalization on the level of 
operators, recently proposed in a slightly different context. We show that in this case the axiomatic 
arrival time distribution of Kijowski is recovered as a limiting case. In addition, it is shown that 
AUcock's complex potential model is also a limit of the physically motivated fluorescence approach 
and connected to Kijowski's distribution through operator normalization. 

PACS numbers: 03.65.Xp, 42.50.-p 



I. INTRODUCTION 

The quantum mechanical analog of the ar- 
rival time of a particle at a given location is 
physically very interesting, and for wave pack- 
ets which are spreading in space this is a highly 
nontrivial subject. It is a particular case of the 
description of time observables in quantum me- 
chanics, i.e., times as random instants - such 
as the arrival times - or durations, e.g. dwell 
or sojourn times. For recent reviews cf. Refs. 
[1, 2]. Difficulties in the formulation of quan- 
tum arrival times and attempts to overcome 
these were presented e.g. in Refs. [3-28]. In 
particular, the lack of a self-adjoint arrival-time 
operator conjugate to the free Hamiltonian lies 
at the core of these difficulties. 

AUcock [3] modeled a simplified detection 
procedure in the region x > by means of a 
complex absorption potential. Because of re- 
flection, he disregarded strong absorption and 
only considered the weak absorption limit, in 
which the detection takes a long time but all 
particles are eventually detected. Under the 
assumption that the measured arrival-time dis- 
tribution was a convolution of an ideal distri- 
bution and an apparatus function he suggested 
for the unknown ideal distribution an approx- 
imate positive expression, obtaining the (not 
semidefinite positive) current density as the ex- 
act solution. Somewhat pessimistically he ar- 
gued that a fully satisfactory, apparatus inde- 
pendent, arrival-time distribution could not be 
defined. 

In contrast, Kijowski [4] (cf. also Ref. [5]) 
pursued an axiomatic approach modeled on the 



classical case and obtained as arrival-time dis- 
tribution at a; = for a free particle of mass m 
coming in from the left with initial state ?A(fc) 
(fc is the wavenumber) an expression which, in 
the one-dimensional case, is given by 



27rTO 



dk i>{k)Vk 



-ihk^t/2m 



(1) 

Surprisingly, this coincides with the approxi- 
mate expression suggested by Allcock [3] . 

Much more recently, the distribution Wk 
has been related to the positive operator val- 
ued measure (POVM) generated by the eigen- 
states of the Aharonov-Bohm (maximally sym- 
metric) time-of- arrival operator [1, 16, 29], 
and this method emphasizes the fact that self- 
adjointness is not necessary to generate quan- 
tum probability distributions. The distribution 
has also been generalized for the case where 
the particle is affected by interaction potentials 
[17-19] and for multi-particle systems [20]. 

Yet, the status of Kijowski's distribution has 
remained unclear and controversial [1, 25, 26]. 
As an ideal distribution, some of its proper- 
ties or of its generalizations have been ques- 
tioned [25] or considered to be puzzling [1], and 
its "operational" interpretation, apart from the 
approximate connection found by Allcock, has 
remained elusive [27]. 

In two recent papers [27, 28], a procedure 
to determine arrival times of quantum me- 
chanical particles has been discussed, which is 
based on the detection of fluorescence photons 
emitted when a two-level atom enters a laser- 
illuminated region. In general, due to partial 
reflection of the atoms by the laser field, not 
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all of them emit photons and hence some go 
undetected. Therefore the measured distribu- 
tion of arrival times is not normalized to one. 
To normalize the distribution, division by its 
time integral was considered ('normalization on 
the level of expectation values'). In some cases 
this gave good agreement with the axiomati- 
cally proposed distribution of Kijowski [4] , and 
parameter regimes where this agreement could 
be found were described. Analogously to All- 
cock 's absorption model, the current density 
could be obtained exactly in the weak laser 
driving limit by deconvolution, and strong driv- 
ing was problematic because of the atomic re- 
flection. The coincidence between the results of 
the simplified complex potential model and the 
more realistic and detailed laser-atom model is 
not accidental and will be explained below. 

Also recently, Brunetti and Fredenhagen [30] 
have proposed a general construction of an 
observable measuring the 'time of occurrence' 
of some event. This construction involved a 
unitary time development and a normalization 
procedure on the level of operators, not on the 
level of expectation values. For this purpose 
they constructed a positive operator on the or- 
thogonal complement of the states on which the 
time of occurrence is infinite or zero and used 
its square root for normalization. This normal- 
ization procedure was in particular applied to 
sojourn or dwell times. 

In this paper it will be shown that normal- 
ization on the level of operators can also be 
applied to the approach to arrival times of 
Ref. [27] which uses spontaneous photon emis- 
sions and, as a technical device, a 'conditional' 
non-unitary time-development. As a result we 
obtain quite simple and explicit expressions. 
In particular, the physically attractive limit of 
strong laser field and fast spontaneous emission 
can be performed and shown to exactly yield 
the axiomatic distribution of Kijowski [4]. 

In the next section we briefly review the re- 
sults of Refs. [27] and [28] and then calculate 
the operator normalized arrival-time distribu- 
tion. In Sections III and IV, fast spontaneous 
emission and strong laser fields arc considered 
in different limits. Finally, a connection be- 
tween the fluorescence approach and complex 
absorption models is exhibited. 



II. OPERATOR-NORMALIZED 
ARRIVAL TIMES 

In the one-dimensional model of Ref. [27], 
a two-level atom wave packet impinges on a 



perpendicular laser beam at resonance with the 
atomic transition. Using the quantum jump ap- 
proach [31] the continuous measurement of the 
fluorescence photons is simulated by a repeated 
projection onto no-photon or one-photon sub- 
space every St, a time interval large enough 
to avoid the Zeno effect, but smaller than any 
other characteristic time. The amplitude for 
the undetected atoms in the interaction picture 
for the internal Hamiltonian obeys, in a time 
scale coarser than St, and using the rotating 
wave and dipole approximations, an effective 
Schrodinger equation governed by the complex 
"conditional" Hamiltonian (the hat is used to 
distinguish momentum and position operators 
from the corresponding c-numbers) 



Hc = ^-ihlm{2\+—e{x)mi\ 



|1>(2|), 
(2) 

where the ground state |1) is in vector- 
component notation , the excited state |2) is 
(°), Q{x) is the step function, and is the Rabi 
frequency, which gives the interaction strength 
with the laser field. 

To obtain the time development under of 
a wave packet incident from the left one first 
solves the stationary equation 

Hc^k = Ek^k, where ^kix) = ,(2), x 

Wk {^) 
(3) 

for scattering states with real energy 

Ek = h^k^/2m 

which are incident from the left (fc > 0). These 
are given by [27] 



^k{x) 



where 
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e''^'" -I- Ri e-^'^'^N 

^ , a; < 

R2 e-"i^ J 

] C+|A+)e^'=+" + C7_|A_)e»'=-^, 

a; > 

(4) 



q = i/fc^ 4- im^/h (5) 
k± = ^k^ - 2mX±/h (6) 

with Img > 0, Imfc± > 0, and where 



A± = (-?:7±iv/72-402)/4 (7) 
1 



|A±) = 



2X±/n 



(8) 



are eigenvalues and eigenvectors of the matrix 

1/0 n \ 

2 In -i-y) 



1 ° The coefficents Ri,R2,C+,C_ fol- 
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low from the matching conditions at x = as 

Rl = [X+{q + k+){k - k^) 

-\^{q + kJ){k-k+)]lD (9) 

i?2 = kik--k+)n/D (10) 

C+ = -2k{q + k-)X-/D (11) 

C_ = 2k{q + k+)X+/D (12) 

with the common denominator 

D = \+{q + k+){k + k-) -X-{q + k-)ik + k+). 

(13) 

By decomposing an initial state as a superpo- 
sition of eigenfunctions, one obtains its condi- 
tional time development. This is easy for an 
initial ground-state wave packet coming in from 
the far left side in the remote past. Indeed, 
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Jo 



dk'ijj{k) *fe(a;) e" 



-ihk^t/2m 



(14) 



describes the conditional time development of 
a state which in the remote past behaves like a 
wave packet in the ground-state coming in from 
the left, with ijj{k), k > 0, the momentum am- 
plitude it would have at t = as a freely mov- 
ing packet. The probability, Nt, of no photon 
detection up to time t is given by [31] 



N{t) = \\^t\ 



(15) 



and the probability density, Il{t), for the first 
photon detection is given by 



m = - 



dN{t) 

dt 



(16) 



For the two-level system under consideration 
one has He — Hi = —i'yh\2){2\, and thus 



(17) 



The integral of the distribution 11 (f) is in gen- 
eral smaller than 1, in fact 



f 

O — ( 



dt n(t) = 1 - iV(oo) 



(18) 



and this was used in Ref. [27] for normalization 

on the level of expectation values. 

In order to employ operator normalization we 
rewrite Eq. (16) in operator form, and to do 
this we go to the interaction picture with re- 
spect to Hq = /2m, 



Hi = e'"°'/^{He- Ho)e-'"'> 



t/h 



^iHot/h^-iH^{t-to)/h^-iHoto/h 



(19) 



where is the conditional time development 
corresponding to H^. Then Eq. (14) can be 
written as 

=e-'^°*/'^C//(i,-oo)|V')|l), (20) 
and Eq. (15) as 

iV(i) = (l|(V;|iV*|V')|l), (21) 

with 

7Vi = C/,^(i,-oo)tC//(t,-oo). (22) 
Similarly, 

n(t) = (i|(V'|fit|V)|i), (23) 

with 

dN{t) 



dt 



(24) 



= 7t^c'(i,-oo)t|2)(2|C//(t,-^). (25) 
In analogy to Eq. (18) we consider the integral 



dtllt = \- No, 



(26) 



and define the operator B on the incoming 
states (with internal ground state) through its 
matrix elements as 

{iMBmi) = {iM\-N^mi). (27) 

The operator 13 can be easily calculated as fol- 
lows. Prom Eq. (14) one sees that for large 

t the second component of ^'(.T,f) is damped 
away and therefore only the reflected wave re- 
mains, 



^{x,t)^ I 
Jo 



dki){k)Ri{ky=' e 



ikx —ihk^t/2m\ 



!)• 

(28) 

for large t. Pulling e"*'*'^ 'Z^™ out from the in- 
tegral as Q-^'Hot/fi- one sees, from Eqs. (20) and 
(22), that 

POO ^ 

[/c(oo,-oo)|V)|l) = / dkip{k)Ri{k)\-k)\l). 
Jo 

(29) 

Taking the scalar product with 
C/c((X), —oo)\(p)\l) one finds from Eq. (27), and 
in k space, 

{l\{k\B\k')\l) = (l - lh(k)Ri{k')y{k - k'). 

(30) 
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Hence, on the incoming states, one can define 
the operator 



^ ^-i/2n 5-1/2 



(31) 



From Eqs. (27) and (24) one sees that 
dttlf^ = 1 and so the probabihty distri- 



bution 



u°-{t)^{iMur\m) (32) 

is normalized to 1. From Eqs. (24) and (30) 
one finally obtains 



/oo /• 
dx / dkdk' i>{k)i){k') (1 - |i?i(A;)n-i/2 (1 - |i?i(A:')|')"'/' 
-CSO J 



X e''^('=^-'=")*/2-<^f (a;)4=^)(a;). 

I 



(33) 



Since |-Ri(/c)| < 1, B is not only a positive op- 
erator but also its inverse square-root exists. 

Operator normalization can be viewed as a 
change in the incident momentum distribution 
4>{k) by a factor of (1 - |i?i(fc)|2)-V2. The ef- 
fect of this factor on a Gaussian wave packet 
is shown in Fig. 1. For mean initial velocities 
of the order of cm/s a single wave packet is 
multiplied by a nearly constant factor. Only 
for very slow particles and $7 ^ 7 a distortion 
of the packet occurs. In this region the am- 
plification of the slow components by operator 
normalization leads to an additional delay of 
n°''(i) compared to U{t). 
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FIG. 1: Operator normalization viewed as change 
of initial momentum distribution. Two Gaus- 
sian momentum wave packets with (t;)i — 2 cm/s, 
{v)2 = 7cm/s, A«i — Av2 = 0.48cm/s, with- 
out (solid hne) and with operator normahzation 
for Q = 0.667 (dashed line) and $1 = 7 (dot- 
dashed line). All figures are for the transition 
6^^3/2 — 6^51/2 of cesium with 7 = 33.3 • 10® s~^. 



III. THE AXIOMATIC ARRIVAL-TIME 
DISTRIBUTION AS A LIMIT 

In this section it will be shown that the 
operator- normalized distribution 11°'' (t) ap- 
proaches Kijowski's axiomatic distribution for 
large 7 and with 7^/0^ = const. We put 
a = -^1 — 40^/7^ and find, for large 7, 



A± 
Q 
k± 
Ri 



R2 




(34) 



to leading order in 7, where the constants Cj 
are given explicitly in the Appendix. Prom this 

one obtains 



(l-|i?i(A:)|2)-i/2(l-|i?i(fe')|2)-i/2 
^ 1 /2m7 



(35) 



and 



7$f (x)*!^^^;) 

hkk' 



277771 



{e(-x)c. 



2g-i(9-9)a; 



^3 I 



-( 



1 + 



2 J I 



'}■ (36) 
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Then n°'^(t) becomes, for large 7 and 0^/7^ = 
const, 



n°~(t) 



2Trm 



J dkdk' iP{k)i^{k')e'^^'''-''"^*^^"'VW 



' —00 
■)2 



(37) 



Inserting q and k± from (34) one sees that the 
expression after x is independent of k and k'. 
One can insert C'i from the Appendix A and 
exphcitly calculate the integral over x, but it 
is easier to note that the term before x is just 
Kijowski's distribution, which is normahzcd to 
1 , and therefore the the expression after x must 
equal 1. 

Thus it follows that 



U°''{t) UK{t) for 7 ^ 00, 'y^/^f = const. 

(38) 

In Fig. 2 it is shown how noN(i) approaches Uk 
for large but finite 7. 




60 80 
t(ns) 



120 140 



FIG. 2: Good agreement of H"'* (circles) with Uk 
(solid line) for large but finite 7, 7 = 107ccsium, 
Q = 0.337. The initial Gaussian wave packet is 
chosen to become minimal when its center arrives 
at a; = (in the absence of the laser) to enhance the 
difference between Uk and the flux (dotted line); 
(v) = 0.9cm/s, Aa; = 0.106 ^im. 



IV. LIMIT OF LARGE O AND 
DECONVOLUTION 

Experimentally, ft is easier to adjust than 7. 
Therefore we also consider the limit of large 
with 7 held fixed. In this case one obtains 



A± 
q 
k± 



^ 2" ~ T 



(39) 



A/fc2 + imj/h, Im g > (40) 



R2 ~ -{l + i)k 



(42) 



mCt' 



to leading order in fl. This yields 

(i-ii?i(fc)|2)-^i-ii?i(A=')n-^ 

and 



(43) 



1 / mfl 



2 V hkk' 
(44) 



/I7 kk' 
2-Km Q. 



{e(-x)2e'(«-«')^ 



+ e(a;)(-ie-'^+^-e-'^- 



X (ie''=+^ -e''=-^)}. (45) 



When integrating over x, only the term 
^-■i{k+-k^)x contributes in leading order of O, 

and this gives 

h 



n°"(i) 



2iTm 

xVkk 



j dkdk' ijj{k)i>{k')e'^ 



7 



7+f(fc2-n 



(46) 



For 7 ^ 00 one again obtains Kijowski's distri- 
bution, but for finite 7 one has a delay in the 
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arrival times. One can try to eliminate this, as 
in Ref. [27], by a deconvolution with the first- 
photon distribution, W{t), of an atom at rest 
in the laser field, making the ansatz 



n°'^(t) = nid(i) * w{t) 



(47) 



for an ideal distribution nid(f). Clearly, W{t) 
has the meaning of an apparatus resolution 
function, similar to Ref. [3] . In terms of Fourier 
transforms one obtains from the ansatz 



(48) 



where [32] 



1 / 7 2\ . 3 , ,2 2 
= 1+ — ^ + - I tuA ^(iuYa ^ 



Prom Eq. (33) one obtains 



7O2' 



^"■^(i,) = J f dx [ dkdk' ^{k)iP{k') (1 - |i?i(fc)n-i/2 (1 _ |i?i(A;')P) 

J —00 J 

x2^s(^'^-^{k'-k'')^W(^)<P^kHx)- 



/N|2n-1/2 



(49) 



(50) 



For large Q one has 1/W{p) ~ 1 + 2^/7. In- 
serting this into Eq. (48) and using Eq. (46) 
yields 



nid(z^) = 



27rm 



/ 



dkdk' il){k)'tl;{k')Vkk' 



x2TrS {u- —(k^ 



fc") ,(51) 



and therefore, for any value of 7 and in the limit 
of strong driving. 



nid(i) = iiKit). 



(52) 



The convergence of Hid to Kijowski's distri- 
bution is shown in Fig. 3. In this example the 
flux, which is a limit of a deconvoluted fluores- 
cence distribution without operator normaliza- 
tion [27], becomes negative. 



V. CONNECTION WITH COMPLEX 
POTENTIALS 

The above approach to arrival times, which 
was based on photon emissions, has another in- 
teresting limit which establishes a connection 
with the complex-potential approach proposed 
by Allcock [3]. Consider now large 7 and O, but 
with n'^/-f = const instead of n^/7^ as before. 
Then a little calculation shows that in Eq. (14) 



the second component V; 



(2) 



7" 



first component goes to 



where 
(t>k{x) 

I! = 

T = 

K = 

Vo = 



dk ^p{k)e~ 



ihk'^t/2n 



'(t>k{x), 



(53) 



e**^^ + Re 

— K 



-ikx 



x<0 
x>0 



k + K 

2k 
k + K 



2imVo 



ImK > 



(54) 



From Eq. (53) one obtains that tA^^^ satisfies 
the one-dimensional Schrodinger equation 

ih^4^^ = (pV2m - iVo&{x))4^^ (55) 



with the complex potential —iVoQ{x). Since 



(2) 



one has, from Eq. (15), 
Nit) = \\4'^\f 



(56) 



and so, from 11 (t) = —dN/dt together with 
Eq. (55), 



V2 while the n(*) = iri (57) 
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FIG. 3: Excellent agreement between the deconvo- 
luted operator-normalized distribution Hid (white 
circles) and Hk (solid line) for large Q = 5OO7. 
Shown is also 11°'^ before deconvolution (dashed 
line). The initial wave packet is a coherent combi- 
nation t/) = 2~^''^(-(/)i + -(/)2) of two Gaussian states 
for the center-of-mass motion of a single cesium 
atom that become separately minimal uncertainty 
packets (with Axi = Ax2 = 0.031 fim, and average 
velocities {v)i — 18.96 cm/s, {v)2 = 5.42 cm/s at 
X = and t — 2 /is) . The flux (dotted) becomes 
negative in some place. 



This is consistent with Eq. (17) since ■y\'tl>t P 
remains finite. 

Eqs. (55) and (57) provide a connection with 
the complex-potential model of AUcock where 
the particle absorption rate is taken as a mea- 
sure for the arrival time. This model is here 
seen to arise as a limiting case from the ap- 
proach of Ref. [27]. It is also obtained by 
considering, somewhat artificially, a position- 
dependent Einstein coefficient, 7(3;) = 76(0;), 
and using an incoming state in the upper level, 
or from the irreversible detector model put for- 
ward by Halliwell [11]. 

The distribution in Eq. (57) is again not nor- 
malized to 1, and it is therefore natural to 
employ an operator normalization. With the 
same arguments as in Section III the operator- 
normalized distribution is obtained as 

n°«(f) = dxj dkdk'jij^)i){k') 

x(i-|i;(fc)|^)-i(i-|i?(fcO|^)-^ 

X T(k)T{k') e»^('=^-'='^)*/2™g-»(«-K')^. 

(58) 

In the limit of strong interaction, Vo ^ 00, one 
again finds that this goes to Kijowski's distri- 
bution, 

n°«(f) ^ n^(i) for Vb 00. (59) 



The advantage of the one-channel model is that 
it provides a simple calculational tool for fur- 
ther, more complicated, arrival time problems 
and that, by simple limits and operator nor- 
malization, it is related to the operational flu- 
orescence approach as well as to the axiomatic 
distribution of Kijowski. 



VI. DISCUSSION 

In Ref. [27] it had been pointed out that from 
the algebraic structure of the arrival time dis- 
tribution in the operational fluorescence model 
it seemed impossible to obtain Kijowski's dis- 
tribution by considering a suitable limit since 
one could not produce the necessary term ^/k. 
This term now arises in the model through an 
operator normalization which corresponds to 
the normalization approach of Ref. [30]. In 
simple, operational terms, this normalization 
can also be viewed as a modification of the 
initial state in such a way that the detection 
losses, due in particular to a strong laser driv- 
ing, are compensated. Our results provide a 
crucial step towards understanding and clarify- 
ing the physical content of Kijowski's distribu- 
tion and, more precisely, establish a set of op- 
erations and limits in which such a distribution 
could exactly be measured. In addition, it has 
been shown in this paper that Allcock's one- 
channel model, which was based on a somewhat 
ad hoc complex absorption potential, is in fact 
a limiting case of the fluorescence model and 
also related to Kijowski's distribution through 
operator-normalization. 

Instead of the operator-normalized expres- 
sion of Eq. (31) one could also consider the 
expectation value of the not manifestly posi- 
tive expression ft/ = i(iJ^^ftt + ftti?~^) whose 
time integral is also 1. Interestingly, in the 
limit 7 00 and 0^/7^ = const this yields for 
the distribution the quantum mechanical flux 
J, discussed in Ref. [27]. 

In this paper we have concentrated on ini- 
tial states with a definite momentum sign, and 
freely moving particles. However, the approach 
can be carried over to a more general setting 
and this will be investigated elsewhere. 



APPENDIX A: EXPLICIT 
EXPRESSIONS FOR Ci(a) 

The constants Ci(a) in Eqs. (34) and (36) 
are given by 



8 



C2 
C3 



2V2a + (1 + a)-V2 _ (l _ a)3/2 



V2aVl - o? + + 1(q - 1) + Vl - q(1 + a) 

\/rT^(x/2 + ^/^^)(a - 1) + ^/^^(^/2 + \/rTS)(a + 1) 
vTTa(V2 + yr^)(a - 1) + \/r^(\/2 + \/rTa)(Q! + 1) 



(Al) 

(A2) 
(A3) 



with a = ./l^^^WPf. 
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